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^ ; RESUMEN 

(S| ; ABSTRACT 

;_^ ' BVR light curves and radial velocities for the double-lined eclipsing binary 

Qh, V1135Her were obtained. The brighter component of V1135Her is a Cepheid 

variable with a pulsation period of 4.22433±0.00026 days. The orbital period 
p>4 ■ of the system is about 39. 99782±0. 00233 days, which is the shortest value 

^N \ among the known Type II Cepheid binaries. The observed B, V, and R mag- 

nitudes were cleaned for the intrinsic variations of the primary star. The re- 
P^ . maining light curves, consisting of eclipses and proximity effects, are obtained. 

C/^ \ Our analyses of the multi-colour light curves and radial velocities led to the 

r^ ' determination of fundamental stellar properties of both components of the in- 

Qh, teresting system V1135 Her. The system consists of two evolved stars, G1-I-K3 

Q \ between giants and supergiants, with masses of Mi=1.461±0.054 M0 and 

M2=0.504±0.040 M0 and radii of Ri=27.1±0.4 Rq and R2=10.4±0.2 Rg. 



< 



U 



^ , The pulsating star is almost filling its corresponding Roche lobe which indi- 

cates the possibility of mass loss or transfer having taken place. We find an 
average distance of d=7500±450 pc using the BVR magnitudes and also the 
^ , V-band extinction. Location in the Galaxy and the distance to the galactic 

^«0 ' plane with an amount of 1300 pc indicate that it probably belongs to the 

thick-disk population. Most of the observed and calculated parameters of 
\Q , the VllSSHer and its location on the color- magnitude and period-luminosity 

^il ' diagrams lead to a classification of an Anomalous Cepheid. 

o' 

fT^ ■ Key Words: stars: binaries: close — binaries: eclipsing — binaries: general 

— binaries: spectroscopic — stars: individual: V1135IIer 

\r2 • 1. INTRODUCTION 

C^ , Double-lined eclipsing binaries are known as astrophysical laboratories. 

Accurate fundamental physical properties of the stars, such as mass, radius 
and effective temperature, are determined directly from their observations. In 
addition they allow us to study internal structure, core or surface convection 
and evolution of the stars. The majority of all kind of s tars occur in binary sys 
tems. Eclipsing bina ry systems include JS cuti stars ( Sovdugan et al. 20061) 



7 Doradus variables (llbanoglu et al.l 120071 ) , RR Lyrae variables (IPrsa et al, 
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2008), Population I Cepheids (jPietrzvnski et alj|2011), Populatio n II Cepheids 
(jWallersteinll2002n and Anomalous Cepheids (jSipahi et alJl2013[ ) which repre- 
sent an important testing basis for theories on stellar structure and evolution 
as well as pulsation mechanisms. Although these variables are located in a 
restricted region in the Hctrzsprung-RusscU (HR) diagram, they have very 
different masses, effective temperatures and chemical abundances. In general 
it is assumed that the components of a binary system are formed simultane- 
ously from the parent cloud. Therefore the components should have the same 
chemical composition at the beginning of their main-sequence evolution. One 
can utilise information from the non-pulsating companion in identifying stellar 
evolution models for pulsating components. 

There arc some metal-poor variables with periods between 0.5 and 3 or 
more days lying above the RRLyrae-type stars which may have entirely dif- 
ferent origin as compared to the normal Type II Cepheids in the globular clus- 
ters. These variables are called as Anomalous Cepheids (hereafter ACs). They 
are found in the general field, globular clusters and nearby dwarf spheroidal 
galaxies. Anomalous Cepheids arc brighter than the Type II Cepheids a give n 
color, and hence they fo llow a different P-L relation ( Soszvhski et al.ll2008bl ). 
Fiorentino et al. I ( 20061) suggested that the ACs are extension of the Type I 
Cepheids to lower metallicities and masses. Moreover, they predict a mass 
range for the ACs as 1.9<M/Mo <3, which may be a product of a mass 
about M<4M0, experienced the helium flash in the past. 

V1135 Her (NSV 10993, 2MASS J18321299-I-1217042, C SC 01032-01378, V =12"'.56, 
B-V=0™. 848) was appeared as a v ariable star in the list of lHofFmeisterl ( 19491 ) 
(see also iKinnunen fc Skifj |2000l )). The hght variabilit y of VI 135 Her has 
bee n attributed to an eclips ing binary of a W UMa type bv lHoffmeisten (|l949l ) 



timated bv lotero fc Wili ( 



and lCoetz fc Wenzel ( 1956). Th e orbital period of the eclipsing pair was es- 



20051 ) as about 40 days and the system classified 



as an EA type with an amplitude of 0.40 mag while it was classified as an 
EW type in the NSV catalogue. Later on it was called as VI 135 Her and 
classified as an EB type eclipsing binary according to the rules of the GCVS 
(IKazarovets et al.ll2008( ). iKhruslovl (|2008l) detected for the first time pulsation 
in the eclipsing binary syst em VI 135 Her based on the NSVS publicly available 
data (jWozniak et al.ll2004l ). He estimated a pulsation period of about 4.2243 
days with an amplitude of ab out 0.2 mag in the R- passband. However, he 
noted that the ASAS-3 data (jPoimanski et al.l 120051 ) did not contradict his 
interpretation. VI 135 Her is an eclipsing binary with a Ce pheid component a n 
analogue of TYC 1031 1262 1, very recently interpreted bv lSipahi et al.l (|2013f ). 
It has the shortest orbital period among the eclipsing binaries having pulsating 
components. 

In this study we present our multi-color photometric and spectroscopic 
observations of V1135 Her. The main aim of this study is to derive the masses 
and radii of the components by analysing the light and radial velocity curves 
using the contemporary methods. Thus, mass and radius of a Cepheid variable 
will be revealed directly from the spectroscopic and photometric observations. 
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TABLE 1 



THE COORDINATES, APPARENT VISUAL MAGNITUDES AND THE 
COLORS OF THE STARS OBSERVED. 



Star 


a 


<5 


V (mag) 


B-V{mag) 


V1135Hcr 
GSC 1032 795 
GSC 1032 1159 


18'' 32"^ 13" 
18 32 12 
18 32 16 


12° 17' 04". 8 
12 11 41 
12 10 34 


12.563 
11.996 
11.730 


0.848 
0.812 
0.935 



We will also discuss the pulsation characteristics and its place in our Galaxy. 



2. OBSERVATIONS 



2.1. Photometric observations 



The photometric observations in the wide-band Johnson BVR system were 
performed with the 35 cm MEADE LX200 GPS telesco pe at Ege Un i versit y 
Observatory. The details of the systems can be found in ISipahi et al.l (|2013l ) . 
The BVR observations in 2008 were obtained on 47 nights between May 30 
and October 22. The observations in 2009 were obtained on 17 nights be- 
tween July 21 and October 8. GSC 1032 795 and GSC 1032 1159 are taken as 
the comparison and check stars, respectively. Some basic parameters of the 
comparison stars are listed in Table 1 . Although the programme and com- 
parison stars are very close in the sky, differential atmospheric extinction 
corrections were applied. The atmospheric extinction coefficients were ob- 
tained from observations of the comparison stars on each night. Moreover, 
the comparison stars were observed with the standard stars in their vicinity 
and reduced differential magnitudes, in the sense variable minus comparison, 
were transformed to the standard system. The standard stars are chosen from 
the lists of Landolt (1983, 1992). Heliocentric corrections were also applied 
to the times of the observations. The V-bandpass magnitude of the system is 
12.56 which is not so bright for a 35 cm telescope. It was observed only three 
or more times in some nights. Therefore we consider that our photometric 
observations were affected from variable atmospheric conditions and position 
of the star on the sky. 

In Fig. 1 we plot the B-, V- and R-passband observations versus the pulsa- 
tion period of the Cepheid variable. As it is seen the dominant light variations 
are originated from the intrinsic variations of the brighter component. The 
below shifted observations correspond to the eclipses. The standard devia- 
tions of each data point are about 0.025, 0.020, and 0.012 mag in B, V and 
R passbands, respectively. The observational data can be obtained from the 
authors. 



2.2. Spectroscopic observations 

Optical spectroscopic observations of VI 135 Her were obtained with the 
Turkish Faint Object Spectrograph Camera (TFOSC) attached to the 1.5 m 
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Fig. 1. The B-, V-, and R-passband light curves for V1135 Her. The ordinates and 
apsis are the standard magnitudes and the pulsation phases. 
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telescope on 13 nights between July 4, 2011 and August 28,2012 under good 
seeing conditions. Further details on the telescope and the spectrograph can 
be found at |http://www. tug. tubitak.gov.tr The wavelength coverage of each 



spectrum was 4000-9000 A in 11 orders, with a resolving power of A/ A A ~7 000 
at 6563 A and an average signal-to-noisc ratio (S/N) was ^120. Wc also 
obtained a high S/N spectrum of the 35 Cyg (F6Ib), 36 Per (F4III), aLyr 
(AOV) and HD 50692 (GOV) for use as templates in derivation of the radial 
velocities. 

The electronic bias was removed from each image and we used the 'crreject' 
option for cosmic ray removal. Thus, the resulting spectra were largely cleaned 
from the cosmic rays. The echelle spectra were extracted and wavelengths 
were calibrated by us ing Fc-Ar lamp source with help of the IRAF echelle 
(JTonrv fc Davislll979f l package. 



The stability of the instrument was checked by cross-correlating the spectra 
of the standard star against each other using the FXCOR task in IRAF. The 
standard deviation of the d ifferences between th e velocities measured using 



FXCOR and the velocities in lNidever et al.l (J2002I ) was about 1.1 kms 



,-1 



3. PULSATION AND ORBITAL PERIODS 
3.1. Pulsation period 



The pulsation period was estimated bv iKhruslovl (|2008l ) as 4.2243 d. Since 
the light variation of the pulsating component dominates in the light curve we 
first attempted to ref ine the pulsation period using the program PERIOD04 



( Lenz fc Bregen 120051 ) . A Fourier power spectrum of all the available V- 



passband data gave a spectral peak at a frequency about /o=0.2367 c/d which 
corresponds to 4.2248 d. This package computes amplitudes and phases of the 
dominant frequency as well as simultaneous multi-frequency sine wave fitting. 
Since the observations were obtained with different instruments in different 
years under dissimilar observing conditions, many systematic observational 
errors and computational errors affected the data. Therefore, we did not at- 
tempt to search for the second or the third frequencies for pulsation if they 
really do exist. As it will be explained in §4 we represented all available V- 
data, subtracting the eclipses, with a truncated Fourier series which includes 
cosine and sine terms up to second order. We then calculated the light vari- 
ation originated from the oscillations of the more luminous component. We 
separated all the data with an interval of about 6 days and determined max- 
imum times by shifting the calculated light curve along the time axis. It 
should be noted here that the shape of the light curve is assumed to be more 
or less constant during the time base of the observations. After the best fit is 
obtained; the times for the mid-maximum light are read off directly from the 
observations and presented in Table 2. Although the observations obtained 
by NSVS and ASAS have relatively large scatters we had to use all the data 
because of limited observations both in time elapsed and the continuous ob- 
servations due to its relatively longer pulsation period. From the time of 
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-960 -840 -720 -600 -480 -360 -240 -120 

Pulsation cycle number 



Fig. 2. The O-C residuals obtained by Eq.l for the pulsating star of V1135 Her 
are plotted against the pulsation cycle number. A linear fit to the data is shown by 
dashed line. 



maximum light and pulsating period given by iKhruslovi (|2008l ) we shifted the 
epoch and using the following elements 



Max{HJD) ^ 2 451 265.430 + 4''. 2243 x E 



(1) 



we obtained the residuals between the observed and calculated times of mid- 
maximum light as well as the number of the elapsed cycles. In Fig. 2 wc plot 
the residuals O-C(I) versus the epoch numbers. The residuals are slightly 
increasing as the pulsation cycle is growing which indicates a correction of 
the pulsation period is needed. A least squares solution gives the following 
ephemeris, 

Max{HJD) = 2 454 653.979(0.148) + 4'*.22433(0.00026) x E (2) 

The standard deviations are given in the parentheses. The pulsation period of 
V1135 Her seems to constant during the observations. T he pulsati o n per iod 
of the Ccpheid is nearly the same with that estimated bv iKhruslovl ( 20081 ). 



3.2. Orbital period 



The orbital period of VI 135 Her was estimated bv IKhruslovl (J2008l ) as 
40 days. We subtracted intrinsic variations of the more luminous star from all 
the available data, thus, the remaining light variations were assumed to be 
originated from the eclipses and proximity effects. We revealed light variations 
due to the eclipses and proximity analysis with a preliminary analysis of these 
data. The shape of the primary eclipse was revealed and compared by the 
observations which fall in the ascending and descending branches of the eclipse. 
Comparing the computed light curve with the observations, we obtained 11 
times for mid-eclipse and presented in Table 3. The low precision is due to 
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TABLE 2 

THE TIMES OF MID-MAXIMUM LIGHT FOR V1135 HER. THE O-C(I) 
AND O-C(II) RESIDUALS WERE COMPUTED WITH EQS.l AND 2, 

RESPECTIVELY. 



HJD-2 450 000 


E 


O-C(I) 


O-C(H) 


Filter 


ReJ. 


1278.573 


-799 


0.190 


-0.167 


R 




1320.551 


-789 


-0.454 


-0.432 


R 




1337.368 


-785 


-0.535 


0.513 


R 




1354.519 


-781 


-0.281 


-0.259 


R 




1358.895 


-780 


-0.130 


-0.108 


R 




1375.938 


-776 


0.016 


0.038 


R 




1414.141 


-767 


0.200 


0.222 


R 




1455.060 


-758 


0.099 


1.010 


R 




1459.708 


-756 


-0.699 


-0.678 


R 




1477.006 


-752 


-0.299 


-0.278 


R 




1481.319 


-751 


-0.210 


-0.189 


R 




1615.087 


-720 


0.492 


0.513 


V 


2 


1619.633 


-719 


0.814 


0.834 


V 


2 


2736.735 


-454 


0.589 


0.602 


V 


2 


2740.946 


-453 


0.576 


0.588 


V 


2 


2854.027 


-426 


-0.400 


-0.388 


V 


2 


2895.924 


-416 


-0.745 


-0.734 


V 


2 


4624.506 


-7 


0.097 


0.097 


V 


3 


4636.930 


-4 


-0.151 


-0.151 


V 


3 


4645.557 


-2 


0.027 


0.027 


V 


3 


4653.978 





0.000 


0.000 


V 


3 


4687.781 


8 


0.008 


0.008 


V 


3 


4700.681 


11 


0.235 


0.235 


V 


3 


5071.970 


99 


-0.214 


-0.217 


V 


3 


5076.351 


100 


-0.057 


-0.060 


V 


3 



Ref: (1) NSVS, (2) AS AS, (3) This study 
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TABLE 3 



THE TIMES OF MID-MINIMUM LIGHT FOR V1135 HER. THE O-C(I) 

AND O-C(II) RESIDUALS WERE COMPUTED WITH THE EPHEMERIS 

GIVEN BY EQS. 3, AND 4, RESPECTIVELY. 



HJD-2 450 000 


E 


O-C(I) 


O-C(II) 


Filter 


Ref. 


1286.75 


-46 


-0.46 


-0.03 


R 


1 


1326.94 


-45 


-0.26 


0.16 


R 


1 


1366.93 


-44 


-0.26 


0.15 


R 


1 


2726.28 


-10 


-0.51 


-0.42 


V 


2 


3126.67 





0.00 


-0.00 


V 


2 


3486.48 


9 


-0.09 


-0.18 


V 


2 


4366.93 


31 


0.61 


0.32 


V 


2 


4646.76 


38 


0.52 


0.16 


V 


3 


4686.79 


39 


0.56 


0.19 


V 


3 


4726.11 


40 


-0.10 


-0.48 


V 


3 


5046.69 


48 


0.57 


0.12 


V 


3 



Ref: (1) NSVS, (2) ASAS (3) This study 

times of minima being based upon very few observations in most cases. The 
O-C(I) residuals were computed using the ephemeris 

Min{HJD) = 2 453 126.2583 -I- 39''.99844 x E (3) 

where the orbital period was derived using the program PERIOD04. In Fig. 3 
we plot O-C(I) residuals versus the epoch numbers. The observations cover a 
short time interval (about 10 years). Since there are four or more nights obser- 
vations on decreasing and increasing parts of the primary minimum obtained 
on JD 2 451286, 2 452726, 2 454646 and 2 454686 a weight of 5 was taken. A 
weighted linear least squares fit to the data gives the following ephemeris, 

Min{HJD) = 2 453 126.678(0.083) + 39'^.99782(0.00233) x E (4) 
The orbital period is very close to that estimated bv iKhruslovl ( 2008[ ). 



4. ANALYSIS 

4.1. Effective temperature of the primary star 

We have used our spectra to reveal the spectral type of the primary com- 
ponent of V1135Her. For this purpose we have degraded the spectral res- 
olution from 7000 to 3 000 by convolving them with a Gaussian kernel of 
the appropriate width, and we have also measured the equivalent widths 
{EW) of the photospheric absorpti on lines for the s p ectra l classification. 



We have followed the procedures of iHernandez et al. I (|2004l ). choosing he- 



lium lines in the blue-wavelength region, where the contribution of the sec- 
ondary component to the observed spectrum is almost negligible. From sev- 
eral spectra we measured i5^WHci-i-FciA4922 = 0.81 ± 0.07 A, i*^WHoi+FciA4i44 = 
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Fig. 3. The O-C(I) residuals obtained by Eq.3 for eclipsing binary V1135Her and a 
linear least squares fit to the data. 



0.33 ± 0.05 A, ^W^caiA5589 = 0.35 ± 0.07 A, EWc^ix6W2 = 0.35 ± 0.07 A, and 
EWr,¥i(n-hi,nd\\4 ^()ri = 1-05 ± 0.02 A. Fr om the calibration relations EW- 
Spectral-type of [Hernandez et al. I (J2004l ) , we have derived a spectral type 
of Gl H-HI for the more luminous star with an uncertainty of about 1 spec- 
tral subclass. In Fig. 4 we compare spectrum of the variable, obtained on 
JD 24 56136 with the spectra of some standard stars. 

We have also observed the variable and comparison stars with the stan- 
dard stars at the same nights. The average standard magnitudes and colors 
of the variable are obtained as <V> =12.563±0.020, <B-V>=0.848±0.040 
and <V-R>=0.551±0.040 mag. Co mparing the location of t he variable in 
the (B-V)-(V-R) diagram given by IPrilling fc Landoltl ( 20001 ) we estimate 
a spectral type of Gl-H-HI which is in a good agreement with that de- 
rived from spectroscopy. The observed infrared colors of J-H=0.316±0.029 
and H-K— 0.121z fc0.028 are obtaine d using the JHK magnitudes given in the 
2MASS catal og JCutri et al. IboOSh. Th ese colors correspond to a G0±2 su- 
pergiant star IPrilling fc Landoltl ( 2000f ) which is consistent with that esti- 
mated both from the spectra and B VR photometry. The effectiye teni pera- 
ture deduced from t he ca. librations of Ide Jager fc NieuwenhuiizenI (|l987l ) and 

Drilling fc Landoltl |2000l ) is 5 280±120K. The standard deviatio ns are esti- 

mated from the spectral-type uncertainty. From the tables given bv lDrilling fc Landolt 
(|2000f ) the intrinsic (B-V) color of 0.82 mag and bolometric correction of 



0.20 mag were interpolated. 
V)=0.028mag is estimated. 



Therefore, an interstellar reddening of E(B- 



4.2. Radial velocities 

To derive the radial velocities of the components, the 13 TFOSC spectra 
of the eclipsing binary were cross-correlated against the spectrum of standard 
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Fig. 4. Comparison of the spectrum of V1135Her with some standard stars with 
similar spectral type but different luminosity. 
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TABLE 4 



HELIOCENTRIC RADIAL VELOCITIES OF VI 135 HER. THE 

COLUMNS GIVE THE HELIOCENTRIC JULIAN DATE, THE ORBITAL 

PHASE (ACCORDING TO THE EPHEMERIS IN EQ. 3), THE RADIAL 

VELOCITIES OF THE TWO COMPONENTS WITH THE 

CORRESPONDING STANDARD DEVIATIONS. 



HJD 2400000+ 


Phase 


Star 


1 


Star 


2 






Vp 


a 


Vs 


a 


55747.3611 


0.5384 


-1.2 


4.3 


-19.9 


5.1 


55751.3423 


0.6379 


8.4 


3.2 


-44.4 


4.3 


55796.3308 


0.7630 


15.0 


0.8 


-61.9 


4.1 


55800.4144 


0.8651 


11.7 


3.1 


-49.1 


3.6 


55835.3383 


0.7385 


14.3 


1.1 


-58.9 


4.4 


55864.2467 


0.4614 


-9.1 


5.2 


11.2 


6.6 


56131.3862 


0.1418 


-17.1 


1.5 


40.1 


4.3 


56133.4240 


0.1928 


-19.4 


1.7 


44.5 


3.4 


56134.3924 


0.2170 


-22.2 


1.3 


49.2 


3.2 


56135.4096 


0.2424 


-18.9 


1.8 


50.9 


3.3 


56135.4505 


0.2684 


-19.9 


1.8 


51.1 


3.8 


56137.3887 


0.2919 


-21.1 


2.0 


49.1 


3.4 


56168.2987 


0.0649 


-10.1 


6.5 


17.8 


7.7 



stars, order-by-ordcr basis using the FXCOR package in IRAF. The majority of 
the spectra showed two distinct cross-correlation peaks in the quadrature, one 
for each component of the binary. Thus, both peaks were fitted independently 
in the quadrature with a Gaussian profile to measure the velocity and errors 
of the individual components. If the two peaks appear blended, a double 
Gaussian was applied to the combined profile using the de-blend function in 
the task. Here we used as weights the inverse of the variance of the radial 
velocity measurements in each order, as reported by FXCOR. 

The heliocentric radial velocities for the primary (Vp) and the secondary 
(Vs) components are listed in Table 4, along with the dates of observations 
and the corresponding orbital phases computed with the new ephemeris given 
in previous section. The uncertainties of the measured velocities are derived 
from different orders. Of course, the velocities of the primary star are slightly 
affected from its pulsation. However, we could not reveal the velocity contri- 
bution originated from the pulsation. Our resolution is inadequate to do this. 
We used HD 50692, Sp.GOV, as a primary standard in the cross-correlation 
which seems to adequate for deriving radial velocities of both components. 
The radial velocities are plotted against the orbital phase in Fig. 5. The 
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0.4 0.6 

Phase 



Fig. 5. Radial velocities folded on a period of 39.99782 days and tiie model. Dots 
(primary) and squares (secondary) with error bars show the radial velocity measure- 
ments for the components of the system 



weight Wi = l/af has been given to each measurement. The standard errors 
of the weighted means have been calculated on the basis of the errors (ai) in 
the velocity values for each order according to the usual formula. The ai val- 
ues are_coni£utedbyM£COR according to the fitted peak height, as described 
by iTonrv fc David (|l979l ) . Radial velocity changes of the Cepheid caused by 
the pulsation are ignored. We did not attempt to decomposition the radial 
velocity measurements of the primary star into the pulsation radial velocity 
and the orbital radial velocity. 

First we analysed the radial velocities for the initial orbital parameters. We 
used the orbital period held fixed and computed the eccentricity of the orbit, 
systemic velocity and semi- amplitudes of the radial velocities. The results of 
the analysis are as follows: e=:0.001±0.001, i.e. formally consistent with a 
circular orbit, Vy^ -4.64±0.42 kms''^, i^i=19.27±0.99 and ii'2=55.88±0.57 
kms~^ . Using these values we estimate the projected orbital semi- major axis 
and mass ratio as: asim=59.39±0.91 R© and q = ^=0.345±0.027. 



4.3. Intrinsic variations of the primary star 

The eclipsing binaries provide critical information about the orbital pa- 
rameters such as orbital inclination, fractional radii, luminosities, ratio of 
effective temperatures etc. If the eclipsing binary is a double-lined binary, 
the masses and radii of the component stars can be determined in solar units. 
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TABLE 5 



THE OBSERVED MAGNITUDES, COLOURS, AMPLITUDES AND RISE 
TIMES FOR THE CEPHEID VARIABLE. 



Passband 


max rain rain 


— max 


M - m 


B 


13.181 13.464 


0.283 


0.388 


V 


12.371 12.600 


0.229 


0.378 


R 


11.840 12.030 


0.190 


0.400 


B-V 


0.810 0.864 


0.054 




V-R 


0.531 0.570 


0.039 





TABLE 6 

FOURIER COEFFICIENTS OF THE OSCILLATION LIGHT CURVES 

FOR VI 135 HER. 



Parameters B V R 

'M, 1.0156±0.0017 1.0081±0.0012 1.0092±0.0011 

Ai 0.1103±0.0024 0.0889±0.0018 0.0715±0.0016 

A2 0.0137±0.0023 0.0132±0.0017 0.0111±0.0015 

Bi -0.0065±0.0024 0.0005±0.0017 0.0050±0.0016 

B2 -0.0175±0.0024 -0.0112±0.0017 -0.0110±0.0016 



Using the inverse-square law one can accurately determine distance to the 
system which is independent of all distance methods. The observed light vari- 
ations are composed of intrinsic light variations of the more luminous star and 
mutual eclipses. First we subtracted all the observations within the eclipses. 
The remaining observations arc phased with respect to the pulsation period 
of 4.22433 days. The light curves of the Cepheid show slow decline and rapid 
rise, i.e., asymmetric light curve with a sharp maximum. The observed mag- 
nitudes, colours and rise times are given in Table 5. The median magnitudes 
are taken as 13.347, 12.514 and 11.964 mag for B, V and R passbands, re- 
spectively. The observed magnitudes are transformed to the flux using the 
median magnitudes. Then, we represented the intrinsic light variations of the 
primary star with a truncated Fourier series. Trial-and-error method showed 
that the observed light curves can well be represented by the second-order 
Fourier series. A preliminary light curve analysis is obtained and the light 
variations out-of-cclipscs, originated from the proximity effects, are revealed. 
These variations are subtracted from the original observations and light vari- 
ations due to the pulsation are obtained. The truncated Fourier coefficients 
are given in Table 6 and the fits are compared with the observations in Fig. 6. 



4.4. Analysis of the light curves 

The intrinsic light variations of the primary star were computed for each 
oscillation phase using the coefficients given in Table 6. After subtraction of 
the Cepheid light changes from the observations we obtained light variations, 
consisting of only from the eclipses and proximity effects. In Fig. 7 the eclips- 
ing light curves in the B-, V- and R-passband are plotted versus the orbital 
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Fig. 6. The B, V, and R passbands light curves of the pulsating primary star and 
their Fourier representations (solid lines) with the coefficients given in Table 6. Note 
that the ordinates are normalized intensities. 



VI 135 HER 15 

phases calculated by the ephemeris given in Eq. 4. The light curve of the sys- 
tem with curved maxima resembles to those /3 Lyrae-type eclipsing binaries. 

We have used the most recent version of the e clipsing binary light curve 



modeling algorithm of IWilson fc Devinnev I (11971) (with updates, hereafter 



W-D), as implemented in the PHOEBE code of iPrsa fc Zwitterl (J2005l) . It uses 
the computed gravitational potential of each component to calculate the sur- 
face gravities and effective temperatures. The radiative characteristics of the 
stellar disks are determined using the theoretical Kurucz atmosphere mod- 
els. The code needs some input parameters, which depend upon the physical 
properties of the component stars. The BVR photometric observations for 
the system were analysed simultaneously. We fixed some parameters whose 
values were estimated from spectra, such as effective temperature of the hot- 
ter component and mass-ratio of the system, which are the key parameters 
in the W-D code. The effective temperature of the primary star has already 
been derived from various spectral type-effective temperature calibrations as 
5 280 K and the mass-ratio from the semi-amplitudes of the radial velocity 
curves as 0.345. A preliminary estimate for the effective temperature of the 
cooler component is made using the depths of the eclipses in the BVR pass- 
bands. Therefore, initial linear and bolometric limb- darkening c oefhcients for 
the pr imary and secondary components were interpolated from Ivan Hammd 



(jl993[ ). taking into account the effective temperatures and th e wave l engths 



of the observations. The bolometric albedos were adopted from iLucv I (|1967( ) 
as 0.5, typical for a fully convective stellar envelopes. The gravity-darkening 
exponents are assumed to be 0.32 for both components, because the stars are 
cool and are assumed to have convective envelopes. The rotational veloci- 
ties of the components are assumed to be synchronous with the orbital one 
and zero eccentricity. At first, we used Mode 2 of the W-D code which is 
for detached binaries with no constraints on the potentials. In this mode the 
fractional luminosity of the secondary component is computed from the other 
parameters via black body or stellar atmosphere assumption. 

The adjustable parameters in the differential correction calculation are the 
orbital inclination, the dimensionless surface potentials, the effective temper- 
ature of secondary, and the monochromatic luminosity of the hotter star. Our 
final results are listed in Table 7 and the computed light curves (continuous 
line) are compared with the observations in Fig. 7. In Table 7 the mean radii 
of the components arc given. The Wilson-Devinney code calculates four di- 
mensions of each component. We give only mean values of the radii which 
are used absolute values of the stellar radii. The distortions on the shape for 
the primary and secondary are Rpole/Rpoint=0.852 and 0.967, respectively. 
The shape distortion of more massive star is larger than that of the secondary. 
The uncertainties assigned to the adjusted parameters are the internal errors 
provided directly by the Wilson-Devinney code. In the last line of Table 7 the 
sum of the squares of residuals (x^), is also given. No reliable solutions could 
be obtained for Mode 4 for the semi-detached binaries with more massive star 
filling its Roche lobe and Mode 6 in which both stars fill their limiting lobes. 
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Fig. 7. The B, V and R passbands light curves, originated only from the eclipses 
and proximity effects, and the computed light curves. The solid lines are for the 
theoretical curves which are also compared by the ASAS V- and NSVS R-passband 
observations in the bottom two panels. 
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TABLE 7 
RESULTS OF THE BVR LIGHT CURVES ANALYSES FOR VI 135 HER. 



Parameters BVR 

r 74.56±0.18 

Te//i (K) 5280[Fix] 

Te//2 (K) 4250±50 

Oi 2.688±0.012 

na 3.368±0.033 

ri 0.4405±0.0025 

r2 0.1682±0.0027 
,, ^V . 0.973±0.004, 0.960±0.004, 0.950±0.004 

X^ 1.470 



We could not obtained any convergent solution from the ASAS and NSVS 
data. The scatter of the observations arc too large to be obtained a reliable 
orbital solution. In the bottom two panels of Fig. 7 we compare our solution 
with the V- and R-passband observations of ASAS and NSVS, respectively. 

We have also applied a procedure, a successive approximation method to 
prevent scatters in the light curves. After subtraction light variation due 
to pulsation the remaining light curves were analyzed and the preliminary 
orbital eleinents were derived. The light variation originated from the prox- 
imity effects removed from the original observations and light variations of the 
pulsating component were revealed. Then the procedure is repeated again. 
Unfortunately the scatters could not be removed. 

5. RESULTS AND DISCUSSION 

Combining the spectroscopic results along with the photometric solutions, 
listed in Table 7, the absolute masses, radii, luminosities and surface gravities 
of the stars are obtained and presented in Table 8. The mass and radius of 
the pulsating component are determined for the first time directly from radial 
velocities and multi-color light curves. They are unexpectedly large compared 
to the Type II Cepheids and small to the Type I Cepheids. The luminosity and 
absolute bolometric magnitude Mboi of each star were computed from their 
effective temperatures and radii. The bolometric magnitud e and effective tem- 



peratu re for the Sun are taken as 4.74 mag and 5770 K ([Drilling fc Landolt 



(|200Cl[ )). After applying bolometric corrections to the bolometric magnitudes, 
we obtained absolute visual magnitudes of the components. Taking into ac- 
count the light contributions of the stars and total apparent visual magnitude 
we calculated their apparent visual magnitudes. The light contribution of 
the secondary star L2/ {Li+L2)=0.027 , 0.040 and 0.050 are obtained directly 
from the simultaneous B-, V-, and R-bandpass light curve analysis, respec- 
tively. This result indicates that the light contribution of the less massive 
component is very small, indicating its effect on the color at out-of-eclipse is 
almost negligible for the shorter wavelengths; 



For the bolometric corrections given bv lDrilling fc Landolti (120001 ) we cal 



culated a mean distance to V1135 Her as (i=7500±450 pc from BVR and JHK 
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Fig. 8. The position of V1135 Her on the color-magnitude diagram. The circles refer 
to the Type II Cepheids in globular clusters and the stars for the binary Cepheids: 
V1135 Her, TYC 1031 1262 1, ST Pup, AU Peg and TX Del. The lines show the 
borders of instability strip for Type II Cepheids. 



magnitudes. The distance from the galactic plane is calculated as 1300 pc us- 
ing its galactic latitude of 9.78 degrees, at least 18 times larger than the mean 
scale height of the Type I Cepheids. The proper motion of this distant Cepheid 
are given by IZacharias et al.l (|2004l ) with too large errors. Therefore, we did 
not attempt to calculate Galactic velocity components. However, its location 
in the Galaxy indicates that VI 135 Her is probably belong ing to the thi c k-disk 
population. This result confirms the earlier suggestion of iHarris et al.l (J1984I ) 
that many field Type II Cepheids are metal rich and have kinematic properties 
similar to the old-disk or thick-disk population. 

In Fig. 8 we plot known Type II Cephei ds in the colo r- magn it ude diagram 



The b inary Type II Cepheids taken from iHarris et al.l (jl984l ) , iNemec et al 



(|l994l ). and lSipahi et al.l (|2013l) were also plotted on the M„-(B-V) diagram. 
VI 135 Her is located among the brightest Type II Cepheids close to the red 
edge of IS, with a luminosity of about 500 times solar. V1135 Her seems to 
an analogue of the binary Cepheids IX Cas with a pulsation period of 9.2 d, 
TX Del with a period of 6.2 d and AU Peg with a period of 2.4 d, ST Pup with 
a period of 18.5 d, and TYC 1031 1262 1 with a period of 4.15 d. They share 
the limited region of the IS of Type II Cepheids. 
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TABLE 8 
FUNDAMENTAL PARAMETERS OF VI 135 HER. 



Parameter 



VI 135 Her 
Primary Secondary 



Spectral Type 


G1(±1)IITII K3(±1)ILIII 


Mass (Mq) 


1.461±0.054 0.504±0.040 


Radius (R©) 


27.1±0.4 10.4±0.2 


Teff (K) 


5 280±120 4 250±50 


Luminosity (Lq) 


517±50 32±4 


Gravity {cgs) 


54±1 129±8 


{vsin i)caic. (km s" 


-1) 34.3±0.6 13.1±0.3 


a(R0) 


61.62±0.95 


Ky (km s-i) 


-4.64±0.42 


z(°) 


74.56±0.18 


q 


0.345±0.018 


d (pc) 


7500±450 



P-L relations for t he Cep heids are studied by ma. ny authors (jStorm et al 
(J2011al ). IStorm et alJ (|2011b[ ). and lStorm et al.l (|2012l ) and references therein). 
The P-L relation slope and zero points h ave long been dis c ussed whether they 
are depended on the metallicity or not. iFreedman et al.l (|200ll ) suggest that 



the zer o point is depended on the metal abundance in the stars. IStorm et al 



(J2011bl ) find identical Milky Way and LMC P-L relation slopes in the near- 
infrared bands by comparing LMC Cepheid P-L relations to their Milky Way 
counterparts. However, zero points show a slight metallicity effect, in the 
sense that metal-poor Cepheids are fainter than metal-rich Cepheids. The 
largest sample of AC s has been collecte d for the LMC by the OGLE HI survey 



( Udalski et al.l 120081) . ISoszvhski et al.l ( 2008a ) presented in the first part of 



the OGLE-HI Catalog of Variable Stars data for 3361 Type I Cepheids in the 
LMC. Of which 1848 are fundamental-mode (F), 1228 first overtone (FO), 
14 second overtone (20), 61 double- mode (F/FO). In the se cond part of the 
OGLE-HI Catalog of Variable Stars ISoszvhski et al.l (|2008bl ) have presented 
197 Type II Cepheids and 83 ACs in the LMC. The T^-band apparent visual 
magnitudes were transformed to the absol ute visual magnitudes My using 
the distance modulus of 18.53±0.02 given bv lPrada Moroni et al.l (J2012l ). The 
Period-Luminosity relations have been determined from a linear regression to 
the the absolute visual magnitudes and log (P) values for the types of Cepheids 
in the form 



Afv 



alogP + h 



(5) 



where My is the visual absolute magnitude and P is the pulsation period in 
days. The coefBcients are presented in Table 9 with their uncertainties. In 



20 



SIPAHI ET AL. 



TABLE 9 

THE PERIOD-LUMINOSITY RELATIONS FOR THE TYPE I 

CEPHEIDS, TYPE II CEPHEIDS AND ANOMALOUS CEPHEIDS IN 

THE LMC. THE FORMAL UNCERTAINTIES ON THE COEFFICIENTS 

AS WELL AS THE NUMBER OF STARS USED AND CORRELATION 

COEFFICIENTS OF THE LINEAR FITS ARE ALSO GIVEN. 



Absolute 
magnitude 



R' 



N 



(mag dcx 



Mv^(CopI, FO) -3.148±0.036 -1.462±0.013 0.86 1228 

Mv(CepI, F) -2.775±0.032 0.953±0.021 0.80 1831 

Mv(CepII) -2.064±0.086 0.409±0.086 0.74 202 

Mv(AC,FO) -2.912±0.384 -1.066±0.092 0.76 20 

Mv-(AC,F) -2.944±0.236 -0.381±0.03.5 0.73 61 



the last two columns the correlation coefficients R^ and number of the stars 
used N are also given. The Period-Luminosity relations for Cepheid variables 
in the LMC are shown in Fig. 9. The individual FO and F pulsators of ACs 
are also plotted in the figure. The FO pulsators are brihter about 0.7 mag 
than the F pulsators at the fixed pulsation period. They lie between the FO 
and F pulsators of the Classical Cepheids (Type I Cepheids). However, F 
pulsators of ACs are located between the F pulsators of the Type I Cepheids 
and Type II Cepheid s in the P-L diagram in contrary to the suggestion by 
Marconi et all (1200J) and ICaputo et al.l (|20M) . They propose that ACs and 
Type I Ceph eids define a commo n re gion in the P-L plan e, confirming the 
suggestion of lDolphin et al.l (|2003[ ) and lCordier et al.l (|2003() that the ACs are 
natural extension of the Type I Cepheids to lower metal contents and smaller 
mass. The pulsating stars in the binary systems are also plotted in the same 
diagram. All of them seem to be located at the same region covered by the 
F pulsators of the ACs. Absolute visual magnitudes, colors and pulsation 
periods are known for only six Type ll binary Cepheids. Th e binarity of the 
distant cepheid QQ Per is doubtful (jWallerstein et al.ll2008| ). 

Individual masses of ACs could be estimated from the period- mass-amplitude 
(P MA) and period-lumi nosity-color relations (PLC). Using the equations given 
bv lMarconi et al.l (120041 ) for the F pulsators of ACs we estimate a mass of 1.3 
M0 from the PMA relation and 1.8 Mq from PLC relation. However, the 
relation for the FO pulsators gives a mass of 1.4 M0. They seem to be in 
agreement within Scr with that found from our analysis . With a mass of 
about 1.461 M0. luminosity of 517 Lq and radius of 27 R© the pulsating 
primary component of V1135 Her cannot be a Type II Cepheid. If it were 
a Type I Cephei d its mass would b e greater than 3.5 M© according to the 
model s given bv IBouo et al.' ("iooO) for Z=0.02 and even for Z=0.004. Re- 
cently |FiOTeiitino_^^Ioiielli (20121 ) suggested that ACs are metal-poor stars 
with the mean mass of 1.2±0.2 Mq. The mass, luminosity, kinematic prop- 
erties and location of VI 135 Her in the C-M and P-L diagrams fulfil most of 
the properties of the Anomalous Cepheids. However, its pulsation period of 
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Fig. 9. The Period-Luminosity relations for the Cepheids in the LMC. Circles and 
dots represent fundamental and first-overtone Anomalous Cepheids. The star sym- 
bols represent binary Type II Cepheids, including V1135 Her. 
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4.22 days is longer than that of the single AGs. The AGs in the LMG have 
pulsation periods generally shorter than 2.5 days. 

Gomparison of its location in the He rtzsprung- Russ ell diagram (HR) , con- 
structed for the metal-poor cepheids bv lCingoldl ( 1985I) . their Fig. 1, indicates 
that V1135Her is in the instability strip of the Type II Gepheids, about 12 
times more luminous than the RR Lyrae stars. It seems to have the same 
luminosity and effective temperature with the Type II Gepheids having pulsa- 
tion periods of about 10 days, closer to the midd le of the IS. When we compare 
with the solar composition models, for example iBaraffe et al.l (|l998f ). the pri- 
mary star appears to be an evolved 5 Mq star and the secondary is consistent 
with an evolved 3 Mq star. 

AGs are found in many nearby Local Group dwarf galaxies but rare in 
globular clusters. There is a general consensus that AGs are central helium 
burning stars with a mass about 1.2 Mq which is too large with respect to 
Type II Gepheids, too small to Type I Gepheids. The origin of their pro- 
genitors is still unclear. There are two suggestions for their origin: a) th e 
evolution of single, met al-poor star with mass < 2 Mq (JBono et al.l (|l997| ) 



Fiorentino et al. I ( 20061)). and b) the evolution of coalescent bina ry systems of 
metal poor stars |Fiorentino fc Monellil |2012l ). Isills et al.l |2009l) ). If the for- 
mer is true they should be younger than the majority of stars in the globula r 
clusters or dwarf sp heroidals in which they are located (JMateo et al.l (|l995[ ). 
Gaputo et al.l (|l999| )). The evolutionary histories of V1135 Her and other close 
binary systems having AGs are very uncertain. 



6. GONGLUSIONS 

We have obtained BVR light curves and radial velocities of the eclipsing 
binary system VI 135 Her of which the more massive star is a Gepheid. The 
astrophysical parameters, mass, radius, effective temperature and luminosity 
of the component stars were obtained directly by analysing the light curves 
and radial velocities. The primary star was classified as a metal-poor Type II 
Gepheids, an analogue of TYG 1031 1262 1. With a mass of 1.46 Mq and radius 
of 27 Rq its structure and evolution are different from those of the RR Lyrae 
and Type II Gepheids. Location of the more massive star on the G-M and 
P-L diagrams indicates that it may belong to metal-poor AGs. We could not 
determine the galactic velocity components due to very large uncertainties on 
the proper motions given in the literature. The distance from the Sun and 
galactic plane are found to be about 7500 and 1300 pc, respectively. All these 
properties and the asymmetric pulsating light curve lead us to classify the 
Gepheid as an Anomalous Gepheid. However, the pulsation period is longer, 
as in the other Type II Gepheid binaries, than those well-known AGs in the 
LMG. 

The pulsating primary star almost fills its corresponding Roche lobe while 
the less massive star fills only 30 per cent. Therefore, mass loss or transfer 
from the pulsating component to its companion has taken place. Both com- 
ponents are located between giant and supergiant luminosity classes and the 



VI 135 HER 23 

separation is only 1.6 times the sum of their radii. The companion should 
unexpectedly affect both pulsation and evolution of the Cephcid. Most of 
the AGs are faint stars and have orbital periods between 40 and 410 days, 
therefore, measurements of the radial velocities AGs in binary systems are 
very difficult. In addition the velocity variations due to the orbital motion 
may be affected by the pulsation depending on the inclination of the orbits. 
In the case of VI 135 Her the orbital inclination is sufficiently large to sepa- 
rate the orbital velocities of both components. If mass loss or transfer from 
the pulsating star has taken place, the orbital period of the system would be 
changed, which can be revealed by the observations to be made in coming 
years. This will yield us some clues about the further evolution of the AGs in 
binary systems. 
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